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Abstract 

We discuss the motion of colloidal particles relative to a two component fluid consisting of solvent 
and solute. Particle motion can result from (i) net body forces on the particle due to external fields 
such as gravity; (ii) slip velocities on the particle surface due to surface dissipative phenomena. 
The perturbations of the hydrodynamic flow field exhibits characteristic differences in cases (i) and 
(ii) which reflect different patterns of momentum flux corresponding to the existence of net forces, 
force dipoles or force quadrupoles. In the absence of external fields, gradients of concentration or 
pressure do not generate net forces on a colloidal particle. Such gradients can nevertheless induce 
relative motion between particle and fluid. We present a generic description of surface dissipative 
phenomena based on the linear response of surface fluxes driven by conjugate surface forces. In 
this framework we discuss different transport scenarios including self-propulsion via surface slip 
that is induced by active processes on the particle surface. We clarify the nature of force balances 
in such situations. 
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I. INTRODUCTION 



Motion of colloidal particles immersed in a fluid can be driven by gradients of concentra- 
tion or pressure even in the absence of net body forces on the particle. In such situations, 
particle motion results from relative motion between particle and fluid induced by surface 
slip As a result, the particle moves relative to the fluid without net force. 

Physical mechanisms that underlie colloidal transport in a fluid can also be used to drive 
the swimming of self-propelling particles. Recently, several scenarios have been discussed by 
which colloidal objects can self-propel. In the case of self-electrophoresis, swimming in a fluid 
is driven by self-generated electric dipole fields, acting at the fluid particle interface 
Phoretic swimmers are driven by a chemical reaction that is catalyzed in an asymmetric 
manner on the particle surface 4, ISl, l6| . The surface reaction generates a concentration field 
near the particle which is asymmetric. The propulsion results from surface slip generated 
by a local concentration gradient at the particle surface. Active processes on a surface also 
propel many cells and microorganisms. An important example is motion driven by motile 
cilia and other active cellular processes [t^. In a coarse-grained description, the beating 
motion of many cilia on a surface effectively generates surface slip velocities that drive 
propulsion j^, 9]. 

Here, we discuss a general theoretical framework to describe the motion of colloidal 
particles in a systematic and controlled way. Our hydrodynamic description is based on 
conservation laws and linear irreversible thermodynamics in a two-component fluid. Exten- 
sion to multicornponent fluids is straightforward. The basic concepts discussed here are well 
estabhshed fllnl, 12|. However, some confusion in the literature calls for clarification [l3| . 

In section II, we review the general hydrodynamic equations of a two-component fluid. 
We express the conservation laws for mass, energy and momentum and identify the conjugate 
fluxes and forces. The hydrodynamic equations follow from a description of dissipative fluxes 
driven by conjugate thermodynamic forces. In section HI, we discuss conditions for which 
concentration gradients and pressure gradients can be generated and the time scales during 
which they persist. A generic theory of dissipative processes associated with slip velocities 
on a solid surface that is in contact with a two-component fluid is discussed in section IV. 
Section V describes the motion of a colloidal particle in pressure gradients and concentration 
gradients. Self-propulsion of a colloidal object due to self-generated concentration gradients 
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as well as due to active surface processes are discussed in section VI. We conclude our work 
with a discussion. 



II. HYDRODYNAMICS OF A TWO COMPONENT FLUID 

In order to define all basic concepts clearly, we first review the hydrodynamics of a two- 
component fluid characterized by a hydrodynamic flow field and diffusive fluxes. Using the 
systematic formulation of hydrodynamic e qua tions based on irreversible thermodynamics 



10|, we follow the discussion presented in 



14j ]. The fluid consists of solvent a and solute 



h with concentrations (numbers of molecules per unit volume) and n^,. The molecular 
masses of the two components are nia and m?,, respectively. We consider the incompressible 
case where the molecular volumes Va and Vb are constant. 

Conservation laws and entropy production. The molecular concentrations satisfy 
the conservation laws 

|i + V.J,^0 (1) 

with i = a,b, where Jj denote particle currents. The mass density of the fluid is p = 
{l — (j))'ma/va + 4>mf,/vb- Here = nf,Vb is the volume fraction of the solute. Incompressibility 
implies that riaVa + UbVi, = 1. Mass conservation can be expressed as 

|^ + V-(pv) = (2) 

where the hydrodynamic flow velocity v is the velocity of the center of mass of local volume 
elements. The particle currents can be decomposed in a center of mass flux and a relative 
flux j = jb = -ja: 

J. = njV + — . (3) 
Momentum conservation is described by the balance equation for the momentum density pv 

dtipVa) - dpa^p = -pgSaz (4) 

where the stress tensor is the (negative) momentum flux tensor and the momentum 
source corresponds to gravitational forces with gravitational acceleration g. The z axis is 
oriented along the vertical direction, opposite to the gravitational field. 
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Dissipation in the system is related to entropy production. The balance of the entropy 
density s reads 

|+V.J, = » (5) 

where is the entropy flux and ^ > is the local rate of entropy production per unit 
volume. Since energy is conserved, the energy density u obeys 

^ + V-J„ = (6) 

where J„ denotes the energy flux. The balance of the free energy density f = u — Ts 
therefore reads 

^ + V-(J„-TJ,) = -r^ (7) 

For an isothermal system (which we consider here) with free energy density /(n^, n^, v, z) = 
{1/2) pv^ + pgz + fo{na, Ub), the total entropy production rate S = J d^rO is given by 

TS = - J (|[^(pv2) + pgz] + ^/^.) + /^^ dA n • (TJ, - J„) (8) 
= j dh{aipu^f)-yyjx) 

where pi — dfo/drii are the chemical potentials of the components. The conservation laws 
and the Gibbs-Duhem relation dP — Uadpa + ribdpf, have been used to obtain the last line. 
Here, denotes the boundary surface of the volume and n a vector normal to the boundary 
pointing out of the volume. The hydrostatic pressure is denoted P. 

The conjugate thermodynamic fluxes and forces are thus the pairs cr^^, Ua/3 and j,— V/i, 
where Ua/s — {daVp + dpVa)/2 is the tensor of velocity gradients and the dissipative part of 
the stress tensor is given by o"^^ = Uap — PSa/s — pVaVp- The relevant chemical potential is 
the difference p — Pb/irif, — Pa/iTia- Expressions for pa and pb in a simple model are given 
in Appendix A. 

To linear order, the dissipative fluxes cr^^, ja depend on the thermodynamic forces Uai3, 
daP as 



ja = -idaP (9) 
1 

3' 



= '^ViUap - l^U^jSap) + rjUjjSap (10) 



Here, the viscosities rj, r] and the dissipative coefficient 7 have been introduced. 
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Hydrodynamic equations. Using Eq. (Q, the solute current can be expressed as 

Jb = -DVub - ^VP + vrib (11) 

Here, D = {'-/ /mh)dfi/dnh\p is the diffusion coefficient of the solute and 7 = 
{•y /mh)dfi/dP\n^ ~ {■y/mh){va/ma — Vh/mb) describes the effects of pressure gradients on 
solute molecules. In the limit of small solute concentration n^Vb <^ 1, the dissipative coeffi- 
cient scales as 7 ~ fimlnb, where is a mobility of solvent molecules. This scaling implies 
that solvent molecules contribute independently to dissipation, no interactions occur. In 
this limit, Eq. (ITTl) becomes a drift-diffusion current for independent solute molecules, see 
Appendix B. 

From the force balance (jlj) and the relation ( fTOl) follows the hydrodynamic equation of 
the barycentric flow field 

r]AY+{f]-'^)V{V■^r) = -VP + pge, (12) 

where is a unit vector in z-direction and inertial forces have been neglected. Incompress- 
ibility of the fluid implies that the molecular volumes Va and Vb are constant parameters. 
Since the density p = (1 — (t>)ma/va + (pnib/vb depends on solute volume fraction 0, the 
divergence 

does not vanish. Here, Ap = rrib/vb — rria/va is the density difference of pure solute and 
solvent. In the incompressible system, the pressure profile P plays the role of a Lagrange 
multiplier function that is determined such that the corresponding flow satisfies the incom- 
pressibility condition ( fT3i) . 

Thermodynamic equilibrium. 

The two-component fluid settles and eventually reaches an equilibrium state. At equi- 
librium, both the flow velocity v and current 3b vanish. The hydrostatic pressure obeys 
dzP = —gp{z). The solute height profile nb{z) at equilibrium satisfies the relation 

p{nb{z), P{z)) = const. (14) 

and is independent on the dissipative coefficients. In the limit of small solute concentration 
n-bVb ^ 1, the solute height profile becomes a barometric distribution nb{z) = ne~^/^ with 
characteristic length I = ksT / [Apgvb), see Appendix B. 
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III. STEADY STATE GRADIENTS 



We now consider a simple geometry to discuss the main features of steady flows and con- 
centration gradients which can be maintained stationary over long times in a two-component 
fluid but are intrinsically nonequilibrium states. A fluid filled channel of height e along the 
X-axis, and infinite extension in y direction, is connected at both ends to reservoirs with 
different solute concentrations or volume fractions (j) = rihVb, see Fig 1. The left reservoir 
is filled with fluid of solute volume fraction 0i up to a height hi, the right reservoir up to 
a height h2 at volume fraction 02- We consider well-stirred reservoirs such that the solute 
volume fraction is uniform in each reservoir. Under such circumstances, is also constant 
in time for short enough times as discussed below. 

For these conditions, the mass densities in the reservoirs are pi,2 = (1 — <Pi,2)^a/va + 
<Pi,2fnb/vb- The pressures in the reservoirs are given by Pi = pigihi — z) + Pq and P2 = 
P25'(^2 — + Pq, where Pq denotes the outside pressure. For /ii 2 ^ e, the pressure gradient 



FIG. 1: Fluid channel of height e and length L between two reservoirs of width w, filled with 
fluid of volume fractions (pi and 02 and heights hi and /i2. In such a setting, gradients of pressure 
or concentration can be maintained over long times. Such gradients can drive the motion of an 
immersed colloidal particle P. 

in the channel generated by the reservoirs is given by 

d,P - g{p2h2 - Pihi)/L (15) 

We consider in the following two complementary cases. 

Case (A): concentration gradient in the absence of barycentric flows. In this 

case, the pressures on each side are balanced, hipi = h2P2, and therefore, d^P = 0, and no 
flow exists, V = 0. In this case, after an initial relaxation process of duration tr — L'^/D, 
the concentration profile becomes a linear gradient 

d,<j)={<P2-<Pi)/L . (16) 
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The corresponding diffusion current is 



{Jb). = -D^^^-^ (17) 

LVb 

This current is maintained on time scales short compared to the equihbration time between 
reservoirs — hwL/{eD), where w and h are reservoir width and height respectively. This 
time can be made arbitrarily long. 

Case (B): flow without concentration gradient. In this case, the volume fractions 
and mass densities are equal in both reservoirs, (pi = 02 and p = pi = p2- For /ii 7^ /i2, a 
pressure gradient and a corresponding hydrodynamic flow exists v = {v^, 0, 0) with V-v = 0. 
The corresponding solution to Eq f|T2|) is 

i»i(^.(.-e) (18) 

In this situation, the solute concentration rib in the channel is constant. A solute current 
exists, {Jb)x = nbVx{z) + jdxP, which consists of a convective part and a dissipative flux 
relative to the fluid flow. This barycentric flow and solute current are maintained during the 
shortest of the times td — hwL/{eD) and tc — wLt]/ {e^ pg). Note again that these times can 
be arbitrarily large. This linear concentration gradient also applies to good approximation 
for weak flows if D/v^ ^ L. For larger flow velocities, nonlinear concentration flelds occur. 



IV. DISSIPATIVE INTERFACIAL PROCESSES: SURFACE SLIP VELOCITIES 

A colloidal particle that is immersed in a two component fluid as described above can be 
set in motion as a result of nonequilibrium conditions provided by either a pressure gradient 
or a solute gradient. In order to determine the velocity of motion of the bead relative to 
the flow, and the corresponding perturbation of the flow fleld, the hydrodynamic equations 
of the fluid are solved with appropriate boundary conditions applied at the particle surface. 
These boundary conditions can be systematically derived by expressing surface dissipation 
at the fluid particle interface and writing generic Onsager relations for conjugate forces and 
fluxes at the interface. 

We consider a solid substrate in contact with a two-component fluid. We use a coordinate 
system with z denoting the distance from the solid surface along the normal direction. In 
the interface, we average all quantities over the thickness d in which properties differ from 
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bulk properties. This procedure is valid when d <^ R, where i? ^ is a local curvature of the 
interface. The interface dissipation rate reads (see Appendix D) 

TS = J dAi^alvt - jtV.fi' + r^A/x) (19) 

Here, the indicies i,j denote directions parallel to the interface. The superscripts s indicate 
that surface fields are considered, vf = Vi{z = 0) and a^^ = aiz^z = 0). Note that we 
have neglected for simplicity surface viscosity and j- is a surface current with units of mass 
per length and time. To keep our discussion simple, we consider here the case where the 
interface equilibrates rapidly with the fluid /i* = fi{z = 0) and material exchange between 
interface and fluid can be negected Vz{z = 0) = and jz{z = 0) = 0. 

In addition to the two conjugate fluxes and forces already discussed, we have added 
here the term r^Afi which describes an active process on the interface which can propel a 
swimmer. An example from biology would be a large number beating of cilia on the particle 
surface. The active process is driven by a chemical fuel (which in a living cell would be 
ATP) with chemical potential difference between fuel and product A/i. The conjugate flux 

denotes the rate of fuel consumption per unit area of the surface. The active process can 
only contribute to net motion generation if the surface has a vectorial asymmetry which for 
example determines the direction of the active stroke of a cilium along the surface. This 
direction is characterized by a normalized vector pi tangent to the surface. 

The three pairs of conjugate fluxes and forces are related to linear order by 

vt = ual-aVijj,' +PiCAfi (20) 
jt = -aat, - YVift' + p.C'A/j (21) 
= Cp^al + Cp^V^^^' + AA^l , (22) 

where the coefficients satisfy Onsager symmetry relations. We observe that there exists in 
general a finite slip velocity at the surface. The coefficient u can be characterized by the 
related "slip length" b = urj, which is the distance from the surface at which an effective 
no slip boundary condition applies. The coefficient 7* is related to surface diffusion and a 
is a dissipative coefficient which couples surface flows to relative fluxes between solvent and 
solute. The coefficients ( and (' describe the coupling of the active process on the surface 
to the two-component fluid. 
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We now consider the case C = cind (' = 0. We furthermore assume that the shp length 
vanishes, z/ = 0, which imphes that there is no shp in the absence of a chemical potential 
gradient. There remains a surface slip velocity driven by chemical potential gradients and 
pressure gradients, which as in the bulk can be expressed as 

v'^ = -nVifil - K'ViP' (23) 

Here k = {a /mb) {dfi/ dnb)\z=o and k' = {a/mb){dfl/dP)\z=o. For a specific interface model, 
the coefficient k has been expressed as k ^ ksTX^ /r] where the length scale A is related to 
the range of the potential describing interactions between solute molecules and the particle 
surface and is of the order of the interface thickness d. 



V. COLLOIDAL TRANSPORT 

We now consider a spherical colloidal particle of radius a which is subject to the solute 
concentration gradient or pressure gradient. We are interested in the particle velocity Vp = 
Vx + Af in x-direction in the laboratory frame where Ai; is the velocity difference between 
particle and flow. These velocities can be determined by solving the Stokes equation f[T^ 
with the appropriate slip boundary conditions on the particle surface. 

In the following, we consider the case where the particle is far from any walls which 
allows us to ignore the effects of boundaries of the channel in which the particle is placed. 
We ignore sedimentation of the particle in the gravitational field in 2;-direction. We discuss 
two cases (A) and (B) described above. 

(A) Concentration gradient in the absence of barycentric flows. In the presence 
of a constant concentration gradient dxUb = (0i — (j)2)/{Lvb), a slip velocity at the particle 
surface is generated. We consider the case of zero slip length 6 = and no active process 
A/i = 0. In spherical coordinates of the particle r, 6, where 6 measures the angle with 
respect to the gradient direction along the x-axis, and r is the radial distance from the 
particle center, the slip velocity v** = Vs{9)eg, is given by 

Vs{0) = -Ksm.{9)dxnb (24) 

where e^ is a unit vector in 9 direction tangential to the sphere. Using the boundary 
conditions veir = a,9) = Vs{9) and Vrir = a,9) = 0, the fiow field around a spherical 
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particle is given by (see Appendix E) 



ve{r,9) 



Vr{r,9) 




(25) 



where we have used for simphcity the condition V ■ v = which is according to Eq. ( |T3|l 
satisfied if Va/rria = Vb/rrib. The relative velocity between particle and fiuid far from the 
particle is 



There is no body force acting on the particle, = 0, see Appendices E and F. The fiow 
field perturbation decays as ~ for increasing r and thus faster than a Stokeslet which 
decays as ~ 1/r and is the signature of a net force. The decay ~ corresponds to a 
source doublet and implies that a force quadrupole is exerted by the particle on the fiuid. 
Here, no force dipole exists. In more general situations, a force dipole also exists, which 
correspond to a decay of the velocity ~ see Appendix E. Note that a force dipole does 
not contribute to propulsion but dominates the far field. Note also, that there are no forces 
due to osmotic pressures acting on a particle in a concentration gradient, see Appendix C. 

(B) Flows in the absence of concentration gradients: In the presence of a pressure 
gradient dxP, a parabolic fiow profile ( fT8|) is generated. If we can ignore the effects of walls, 
and if no surface slip occurs, the relative velocity of the particle with respect to the fiow 
At>i = Vp — Vx can be determined by Faxens law as (see Appendix E) 



Surface slip induces an additional component Av2 to the relative motion between fiuid and 
particle. For vanishing slip length 6 = 0, the only contribution to the slip is = —n'ViP'^. 
The boundary conditions in a coordinate frame co- moving with the particle are Vr{0,r = 
a) = and vg{6,r = a) = —k' sin{6)dxP- The fiow field perturbation generated by surface 
slip is given by Eq. ([25]) with Av = Av2 = {2k' /?,)dxP. 

Since k' expresses the momentum transfer to the particle within the interfacial layer of 
thickness ci, one expects ^ —(P/rj, where the length d is of order d. By superimposing 
the fiow field (ITSl) with the perturbation by the particle At>i and At>2, the overall particle 



Av = {2K/3)dxnh 



(26) 



Avi 



Sa^dxP 
4r) 



(27) 
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velocity is 



g{p2h2 - pih) 
L 



(28) 



81] 12r] 3?7 



which is entirely driven by the pressure gradient. No net body force acts on the particle, 
see Appendix F. Note that the convective term is large compared to Avi which is in turn 
large compared to Av2- The existence of Avi is yet important since it can lead to separation 
of particles according to their size. 

VI. SELF-PROPULSION 

Propulsion by self-generated concentration gradients. A chemical reaction cat- 
alyzed on the particle surface can generate local concentration gradients which propel the 
particle 0, . To describe situations where the particle catalyzes a reaction involving two 
solvent species, at least a three-component fluid description is required. However, the basic 
physics of self-propulsion can be captured by our two-component fluid if we assume that a 
surface reaction transforms molecules of type a into molecules of type b. In this case, we 
consider two components with equal molecular masses rria = to satisfy mass conserva- 
tion. The molecular volumes Va and can in general differ. A local reaction rate S per unit 
area at which a molecules are transformed in h molecules on the particle surface implies the 
boundary conditions for the molecular fluxes normal to the particle surface 



in the reference frame where the particle is at rest. The corresponding boundary condition 
for the center of mass velocity is Vz{z = 0) = 0. 

The nonlinear convective term in the flux Eq. (fTTl) . couples the concentration and flow 
fields. In the limit of small Peclet number Pe = Av a/D, we can neglect this convective 
nonlinearity. The stationary concentration field is then solution to the diffusion equation 
= 0, 0]. For the simple choice of an asymmetrically distributed reaction rate S{9) = 
So cos{9), the solution for the boundary conditions specified above is given by 



where denotes the concentration far from the particle. This concentration field induces 




(29) 




(30) 



11 



according to Eq. (l23l) the surface slip velocity 

^sW = |§sin(e) (31) 

which is independent of particle radius. The corresponding hydrodynamic flow field for the 
case Va = Vb for which V ■ v = is given by Eq. fl25|) with 



= . (32) 

The pressure P is constant and there is no net body force acting on the particle. 

For finite Peclet number, there exist no simple solutions to this nonlinear problem since 
the convection velocity in Eq. flTTl) is given by the hydrodynamic flow field described by Eq. 
f ll2p . The latter in turn is coupled to the concentration field Ub via the surface slip driven 
by the chemical reaction. Here, we focus on the scaling behavior of the propulsion velocity 
for large Pe. The hydrodynamic flow relative to the particle is of order Av. This defines a 
length scale i ~ (Da/ AvY^"^ <^ a which characterizes the thickness of the boundary layer 
near the particle in which the concentration field is different from n'^ by a variation of order 
Arib — iSo/D. The resulting slip velocity is of the order Av ~ KAub/a. By combining these 
expressions, we obtain 

^2/3 o2/3 

The velocity thus increases for large Pe less than linearly and becomes dependent on the 
particle size a. For Pe = 1, the relations fl32|) and (l33l) match except for a dimensionless 
prefactor. 

Propulsion by active surface processes. Many microorganisms swim by using a 
large number of cilia attached to the surface which generate periodic beating movements 
that are driven by molecular motors that consume a chemical fuel [3]. In a coarse-grained 
picture where the cilia are active elements within the interface between swimmer and fluid, 
the motion of the cilia effectively generates a surface slip velocity of the flow. In our generic 
description this surface slip is captured by the term ^ Pi^Afi in Eq. fl20|) . Here, the 
tangent vector pi describes the direction along which the cilium generates a flow and ( is 
a coupling coefficient between the free energy Afi driving the motors and the surface flow. 
In this scenario, the slip pattern on the surface of the swimmer is determined by the ciliar 
beat direction and strength. 
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For given surface slip pattern and corresponding hydrodynamic flows with V ■ v = 0, the 
velocity of the swimmer can be determined without explicit calculation of the flow field. For 
a spherical particle it is given by 1^ 

Similarly, the rotation rate and axis (which can exist for a surface slip velocity field lacking 
axial symmetry) is described by [l^ 



n = - — - / ciA n X (35) 

Sna-^ J 

where n is a unit vector normal to the surface and f2 points in the direction of the rotation 
axis. Again, this translational and rotational motion occurs relative to the fluid without a 
net body force and torque acting on the swimmer. For the simple example with pi a unit 
tangent vector in 9 direction and C = Co sin(6'), we have Av = {2/3)CoAfj, and again the flow 
field perturbation of Eq. (l25l) . 

VII. DISCUSSION 

We have presented a generic description of colloidal transport to clarify the force balances 
involved and the role of interfacial slip. The simplest form of transport occurs if relative 
external body forces are applied to a particle. External body forces, such as those due 
to gravitation, correspond to source terms in the momentum balance (jlj). As a result, the 
particle moves at a speed Av relative to the fluid and experiences Stokes friction = 
GirrjaAv. The corresponding perturbation of the fluid flow at large distances is given by a 
Stokeslet which decays as 1/r. 

In the absence of external fields, a colloidal particle in a fluid is force free, = 0, even if 
concentration or pressure gradients exist. In this case a colloidal particle nevertheless moves 
relative to the fluid at a velocity Av if there is a slip of the flow at the surface of the particle. 
Such slip is generated by surface dissipative phenomena We have shown that Onsager 
relations on the solid surface determine the boundary conditions for the hydrodynamic 
equations in the bulk. This description can account for a variety of phenomena, including a 
slip length, slip induced by concentration or pressure gradients as well as slip due to active 
processes on the particle surface, within a unified framework. 
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Our generic description can be generalized to electric fields, which requires to replace 
chemical potentials by electrochemical potentials and to include electrostatics in the free 
energy density. Effects in the presence of electric fields include electrophoresis of charged 
particles. Since a charged particle is screened beyond an electric double layer of counter 
ions, it is effectively neutral. Electrophoresis is thus another example where no force acts on 
a particle (including the layer of counter ions) and electrophoretic motion results thus from 
surface slip 

Our arguments are relevant for mechanisms of self-propulsion of colloidal particles. Self 
propulsion implies that motion occurs in the absence of externally applied forces. This is 
possible if a particle self-generates a surface slip velocity by active processes on or near the 
particle surface. These include the action of cilia and flagella in the case of swimming mi- 
croorganisms Ql or the generation of a concentration gradient by surface chemical reactions 
6|. The main effect of the surface slip is to generate a relative motion between particle and 
fluid. In the general case, this relative motion is associated with a flow field perturbation 
that decays as ~ corresponding to a force dipole exerted by the particle on the fluid 
which however does not contribute to propulsion. If no force dipole exists, such as in the 
case described by Eq. f l25|) . the flow perturbation decays as ~ This flow perturbation 

corresponds to a force quadrupole acting on the fluid. In all cases there is no net external 
force acting on the particle. 

In a recent publication it has been claimed that self-generated concentration gradients 
generate a force on a particle by osmotic pressure gradients which is balanced by Stokes 
friction ISj. Clearly, no force is exerted on a particle in a concentration field by osmotic 
pressure, see Appendix C. The scenar io p roposed in [13l] violates momentum conservation, 
see Appendix F. The work of Ref. [l^ furthermore ignores the effect of hydrodynamic 
flow perturbations on the concentration gradients, when those cannot be neglected at finite 
Peclet numbers. Our generic description of colloidal transport can serve to clarify these and 
related points in a systematic way. 

We thank Ramin Golestanian for stimulating discussions. F.J. thanks Friederike Schmidt, 
Holger Stark and Andrej Vilfan for similarly stimulating discussions on the hydrodynamic 
flow perturbations generated by colloidal particles and swimmers. 
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APPENDIX A: ENTROPY OF MIXING AND CHEMICAL POTENTIALS 



The properties of the chemical potentials can be discussed using a simple model for a 
two-component fluid with a free energy in the rest frame Fo{Na, Nf,, V) — Vfo{Na/V, N^/V), 
where V denotes volume, with 

fo{na,nb) = kBT(naln hn^ln ) + ^{naVa + ribVb - 1)^ (Al) 

riaVa + ribVb riaVa + TlbVb 2 

Here, the first term describes the entropy of mixing of two components with molecular 
volumes Va and Vb, the second term describes the compressibility of the fluid by the coefficient 
X- Interactions between the two components are neglected. The chemical potentials /li — 
dfo/dui are given by 

/ N 1 rr ''^o.Va , nb{vb -Va)\ , . 

fJ'airia, rib) = ksT In ■ \ ■ + VaX{naVa + ribVb - 1) (A2) 

V riaVa + UbVb UaVa + UbVb J 

pibiria: rib) = ksT In ■ \ ■ + VbX{naVa + ribVb - 1) (A3) 

\^ naVa + ribVb riaVa + UbVb J 

The pressure P = -(9Fo/9y)|jv„,iVi = -/o + /^a^^a + l^'brib is 

P(n„, Ub) = -{{riaVa + ribVb) - 1) (A4) 

Using the pressure, the solvent density Ha can be eliminated. In the incompressible limit of 
large x, riaVa + ribVb — 1 and we obtain 

jJiairih, P) - kBT{hl{l - HbVb) + nb{Vb - Va)) + Pva (A5) 

Hb{nb,P) - kBT{\n{nbVb)+na{va-Vb)) + Pvb (A6) 



APPENDIX B: LIMIT OF SMALL SOLUTE CONCENTRATION 

The solute flux is driven by gradients of the chemical potential difference Jx = iib/rub — 
i^alrua, 

3b^ -—Vjl + UbV (Bl) 
nib 

In the limit UbVb <^ 1 of small solute concentration, 

fl{nb, P)^^ HubVb) +(^-^)p (B2) 
rUb Krub rUaJ 
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For small n^, solute particles become independent of each other and dissipation takes place 
independently for each solute particle. Therefore, 7 ~ ^mfn;,, where ^ is a mobility per 
solute molecule. We thus find 

Jf, ~ -DVuh - ^n^VP + n^w (B3) 

where 

7 - ^mb{vb/mb - Va/rria) 

D ~ ^keT (B4) 

In the limit of small HbVb, the pressure gradient is approximately constant dzP = —pg — 
—grria/va- For v = 0, the height profile is 

Ub = ne-'/^ (B5) 

with i = Dva/{'yiTLag). Using Eq. flB4l) . we obtain i = kBT/{ApgVb) and flB5l) is the 
barometric height distribution. 

APPENDIX C: OSMOTIC PRESSURE 

Osmotic pressures are a consequence of a semipermeable interface which separates the 
fluid in two compartments (1) and (2). The solvent passes this interface, which is imper- 
meable to the solute. As a consequence, across the interface the chemical potential of the 
solvent is balanced, /i^^-* = However, the solute chemical potentials do not balance 

(1) / (2) 

The solvent chemical potential is according to Eq. flA5l) in the limit of small UbVb given 
by fia — —ksTubVa + Pva- The balance of solvent chemical potentials implies the existence 
of an osmotic pressure difference across the semipermeable membrane 

p(2) _ p(l) ^ ^^2^(^(2) _ ^(D) . (CI) 

Note that the hydrostatic pressure difference appears only after the balance of the chemical 
potential of the solvent is reached and that the corresponding momentum source is provided 
by the membrane. 
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APPENDIX D: FORCE BALANCE AND DISSIPATION AT AN INTERFACE 



We consider dissipation and force balances in an interfacial region of thickness d between 
two phases in which material properties differ from those in the two bulk phases. The 
coordinate normal to the interface is denoted z. A relative slip velocity can occur at an 
interface. The local center of mass velocity tangential to the interface aX z = ±d/2 is 
denoted vf, respectively. Dissipation due to interfacial slip 



vt = vt-vr (Dl) 



/■'^/2 8 V 1 

TSc^ I dA dz^ai.c^- I dA a^v^ (D2) 



can be expressed as 



-d/2 2 2 

The interfacial shear stress crf^ = ea^^ + {l — e)a~ is a weighted average of a^. The value of 
< e < 1 depends on the internal structure of the interface. The conjugate thermodynamic 
variables are thus ff and af^. The corresponding Onsager relation reads 



u{eat + (1 - t)a-) (D3) 



where v is the corresponding dissipative coefficient. 

The force balance in the interfacial region d^aap = implies 

rd/2 

-d/2 

This implies the interfacial force balance 



rd/2 

/ dz + d,<j,,) = (D4) 

J —dl2 



(D5) 



where the interfacial tension (for isotropic stresses in the tangent plane) is 

1 rd/-2 

S = - / dz akk (D6) 

2 J -d/2 

In the absence of interfacial tension gradients SjE, the shear stress is continuous across the 
interface af^ = af^ = a^^. and the slip velocity is simply given by vf = vaf^ as in Eq. (120|) . 
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APPENDIX E: HYDRODYNAMIC FLOW FIELDS WITH AXIAL SYMMETRY 



Solutions to the Stokes Eq. (fT2l) for axisymmetric incompressible flows with V ■ v = 



can be expressed using the stream function ip 
field is related to the stream function ip{r,9) by 

1 dip 
^ sin 9 36 



15| . In spherical coordinates, the velocity 



(El) 



ve = (E2) 
r sm & or 



The stream function satisfies the differential equation E^^ip = 0, where 

(9^ sinO d 1 d 



^2 , ( sin^ d I d \ , 



Simple solutions are given by 

^ = sin^ e{Air^ + A2r^ + Agr + ^) (E4) 

where Ai, . . . , are constant parameters determined by boundary conditions. The corre- 
sponding pressure field is given by P = — cos 0(2Oy4ir + 2A^/r'^) + P°°, where P°° is the 
pressure far from the particle. The body force acting on the particle which is balanced by 
forces exerted by the hydrodynamic flow is = —Sin] A3 [isl. 

In case (A), we determine a solution of the form given by Eq. ( 1E4I) with V0{r = a) = Vs{6) 
and Vr{r = a) = in the reference frame moving with the sphere. For large r, we require 
motion at constant velocity Av in negative x-direction, Vr — —Av cos{6) and ve ~ Av sin 6. 
From the latter conditions, it follows that Ai = 0, and 2A2 = —Av. Because there is no 
external force acting on the particle, f '^ = and thus ^3 = 0. The boundary conditions on 
the particle surface imply A4 = —A2a^ and 3A2 = —{K/vb)dx(f)- The corresponding flow is 
given by Eq. ( 125|) . The perturbation of the flow velocity decaying as ~ corresponds to 
a source doublet 18 1 and implies that a force quadrupole is exerted by the particle on 
the fluid. 

In case (B), we superimpose the flow ( IT8l) driven by an applied pressure gradient with 
the velocity Af 1 due to the perturbation of the Poiseuille flow by the sphere and the relative 
velocity Av2 due to interfacial slip. The velocity Af 1 can be estimated from Faxens theorem 
which expresses the force on the particle as 

= Qnvaiv, - vl) + vra^V^f ° (E5) 
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where denotes the unperturbed parabohc flow fleld. No force acts on the particle, = 0, 
which determines the relative velocity Af i = Vp — v^. This flow is superimposed with a flow 
driven by slip boundary conditions (12^ as described for case (A). 

If a net body force acts on the particle with no slip, the boundary conditions are 
Vr{r = a) = and V0{r = a) = which require 3^4 = A^a? and A2 = —(2/3)^13. From the 
asymptotic behavior, it follows that Ai = and 2A2 = —Av. This implies = GnrjaAv 
and the flow fleld is given by 

ve{r,e) = At; (^1 - - sin(^) (E6) 

In addition to the force monopole this flow contains again a contribution from a source 
doublet implying a force quadrupole. 

The equation E^tp = for the flow also has the solution 

^ = - cos ^ sin^ e{Bi + Es/r^) . (E7) 

The corresponding flow fleld is Vr = (3 cos^ — l)(i?i/r^ + i?2/r^) and vq = B2 sm{26) / {2r^) 
[9]. The condition Vr{r = a) = that the radial flow vanishes on the particle surface imposes 
B2 = —Bio?. The decay of the radial component proportional to ~ corresponds to 
a stokes doublet which implies the action of a force dipole on the fluid 18|. This force 
dipole dominates in the far fleld over the force quadarupole. For an arbitrary distribution 
of surface slip ve{r = a) = Vs{0), a force dipole exists in general. The force dipole vanishes 
by symmetry if the surface slip is a symmetric function Vs{6) = fs(7r — 6) such as is the case 
for Vg ~ sin(6') described by Eq. fl2S]) . Note that the flow perturbation corresponding to a 
force dipole does not contribute to propulsion by symmetry. 



APPENDIX F: GENERAL CONSIDERATIONS CONCERNING THE EXIS- 
TENCE OR NON-EXISTENCE OF A STOKESLET IN A VELOCITY FIELD CAR- 
RYING A PARTICLE 

Consider a particle of any shape, including arbitrary topological genus, immersed in a 
multi-component fluid flowing in a container of complex geometry possibly of non trivial 
topology. The fluid is submitted to an external force of density g^'"^^^ and the particle to the 
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force density g'^'^^*. The momentum flux in the fluid is characterized by a stress tensor Gap- 
The total force acting on the particle reads: 

/r = / gT'^'dV-r [ a^pdAp (Fl) 

Here, the integration volume Vp and surface Sp refer to the particle and the surface elements 
dAp are oriented to point outward from the particle. The particle dynamics is mXa = f^^, 
where Xa is the particle position and m is the particle mass. In the Stokes limit where inertial 
terms can be neglected or in stationary conditions /*°* = 0, which implies that in inertia 
free regimes the total force acting on a particle in a fluid vanishes. Momentum conservation 
in the fluid implies 

gff"^ + dpa^p = (F2) 
where inertial forces have again been neglected. In integral form, this can be expressed as 

/ a^pdAp+ [ a^pdAp = - f g^f^'dV (F3) 

Here the integration is over the particle surface Sp and an arbitrary surface S enclosing a 
volume V which includes the particle in the fluid. The surface elements are oriented to point 
to the outside of the fluid volume V — Vp between the surfaces Sp and S. The (vanishing) 
total force now reads 

/r = / + / a^pdAp + / gff'-'dV (F4) 

t/Vp v S v — Vp 

or equivalently, 

/r = / {g'r' - 9t''')dV + j a^pdAp+ [ gff-'dV (F5) 

t/ Vp v 5^ v V^ 

Consider now the stress (j° in the fluid in the absence of the particle (every other condition 
being kept identical). For any volume and corresponding surface we have 

/ gff-'dV ^ - j a%dAp (F6) 

•/ V^ J s 

Therefore, for any surface S enclosing the particle in the fluid, 

fT = L igr' - 9^f''')dV + [ (a^p - a%)dAp = (F7) 

J Vp J o 

We thus find that the relative external body force 

/S = / {g'r' - 9t''')dV (F8) 

t/ Vp 
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which includes an Archimedian correction to the net body force, is balanced by the stresses 
exerted by the perturbation of the flow field on the particle 

Fa = - li^cp - alp)dAp (F9) 

(for the choice 5* = Sp). 

Eq. (]F9|) is valid for arbitrary S enclosing the particle. Therefore, for colloidal particles 
far from the boundaries of the fluid, the total momentum flux through any closed surface is 
constant and equal to This implies that the stress field perturbation cTq,/? — a^^ ~ 
for large r, where r is the distance from the particle center. Since aa/3 — 0"°^ ~ {daVp + dpVa), 
the barycentric velocity perturbation scales as |v| ~ 1/r which implies the existence of a 
Stokeslet. If the relative external body force vanishes, but a relative external torque acts 
on the particle, a similar argument using angular moment conservation implies that in this 
case the stress scales as ~ and the velocity field like ~ 

If no external field acts on either the particle or the fluid, there is neither external body 
force no torque. The far field of the hydrodynamic flow is dominated by a force dipole or 
Stokes doublet with stress ~ and velocity ~ If the force dipole vanishes, the 

stress scales as ~ and the velocity as ~ As we have illustrated in this work, 

the absence of external forces and torques does not mean that there is no motion between 
particle and fluid. 
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